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Based on our experience in kinetic modelling of coupled multiple metabolic pathways, 
we propose a generic rate equation for the dynamical modelling of metabolic kinetics. 
It is symmetric for forward and backward reactions. It's Michaelis-Menten-King-Altman 
form makes the kinetic parameters (or functions) easy to relate to experimental values 
in database and to use in computation. In addition, such uniform form is ready to ar- 
bitrary number of substrates and products with different stiochiometry. We explicitly 
show how to obtain such rate equation rigorously for three well-known binding mecha- 
nisms. Hence the proposed rate equation is formally exact under the quasi-steady state 
condition. Various features of this generic rate equation are discussed. In particular, for 
irreversible reactions, the product inhibition which directly arise from enzymatic reaction 
is eliminated in a natural way. We also discuss how to include the effects of modifiers 
and cooperativity. 

Keywords: Generic Enzymatic Rate Equation; Metabolic engineering; Systems Biology 
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1. Introduction 



The advances in modern biology require large-scale mathematical modelling 
[TJ |5J 121 H] . One such important task is the kinetic modelling of coupled metabolic 
pathways [51 |SJ |7J ISl IE EHl IHl 112 , to understand how an organism lives. While it is 
true that biological processes are based on well-studied chemical reactions, a mod- 
eler immediately encounters several formidable difficulties. First, kinetic equations 
describing enzymatic reactions are often complicated due to numerous parameters 
|131 114j. even without considering stochastic effects Moreover, each enzymatic 
equation appears different from others depending on the enzyme mechanism in- 
volved. Therefore, modelling a simple biological process with relatively detailed 
chemical and biological information can be a daunting task [T5] [T7] . Thus, those 
rigorous chemical reactions in a coupled metabolic network do not lend themselves 
easily to metabolic simulations with many reactions. Nevertheless, many databases 
such as Brenda |TS1 US], KEGG [2D], and MetaCyc [2T1 US have been document- 
ing our progress in knowledge of enzyme behavior. If we assume that this difficulty 
could be overcome with great care and effort in using the appropriate rate equations 
with the aid of enzyme databases, there exists a second difficulty. Most chemical 
reaction parameters are not measured under the living conditions, how does one 
know they are applicable to a real organism? If not, how could one make the ap- 
propriate and necessary adjustments in a complicated though exact rate equation? 
More seriously, even if those in vitro parameters are relevant, it is not possible to 
measure all of the rate constants in the rigorous rate equation, whose number can 
easily be tens of thousands for multiple coupled pathways [23j [24] in ever changing 
physiological conditions, the third difficulty. The question naturally arises: Is there 
a way to address those difficulties? 

In this paper, based on our kinetic modelling experience we introduce a generic 
rate equation capable of reproducing exactly the full rigorous rate equations irre- 
gardless of the enzyme mechanism for the full rate equation. The equation takes into 
account thermodynamic constraints. It generalizes easily to any reactions containing 
an arbitrary number of substrates and products. The effects of modifiers (activa- 
tors and inhibitors) as well as cooperativity can be included. Thus, it appears to be 
capable of addressing the above difficulties. 

The outline of the rest paper is as follows. In the next section. Sec. |2| we intro- 
duce the generic enzymatic rate equation. In Sec. [3] we show how the full rigorous 
mechanistic rate equation can be written exactly into the form we proposed using 
three important cases as examples. In Sec.|4] we discuss the case when cooperativity 
is involved. In Sec. [5] we discuss a useful ansatz for the functions /i and /2 and also 



thermodynamic considerations in the generic equation, Eq. (2.2 1. We also discuss 
relationship to another proposed rate equation that is different from ours and how 
to include the effects of modifiers. We conclude in Sec. [6| 
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2. Proposed Generic Enzymatic Rate Equation 

In this section we introduce the generic rate equation. 

To further illustrate difficulties and to demonstrate usefulness of an answer to the 
questions posed in the introduction, we consider the well known Michaelis-Menten 



equation (see Eq. (3.11 1). Once we are dealing with bisubstrate enzymes beyond the 
simple Michaelis-Menten equation [2S] for a single substrate and product, the rate 
equations becomes unwieldy due to the rapid growth in the number of parameters 
in the equation. This is despite the fact that the Michaelis-Menten equation has 
been shown to hold even for a fluctuating single enzyme ^26j . Is it possible to fully 
measure the 18 terms in the steady state |571[5H] rate equation for a random Bi-Uni 



mechanism (c./. Eq. (3.13 1)? The answer has been no so far experimentally. It has 
been recognized that enzymes have not been characterized to such great details in 
biochemical literature [35]. Various rate equations in bisubstrate enzyme reactions 
alone can be found in |30| . However, as pointed out in Ref. [31 , the primary concern 
with a rate equation used in a metabolic simulation isn't to account for the fine 
difference in kinetic behavior predicted by different mechanism. Instead, we need an 
equation that is amenable to a wide variety of enzymatic reactions but is still capable 
of describing the response to substrate variations and product concentrations |29j . 



To this end, a generic rate equation (c./. Eq. (2.2|) will be very useful since it 
can represent exactly the original mechanistic rate equation. Irregardless of whether 
the mechanistic details of a reaction is known, such a generic rate equation can 
be adapted to all the reactions that we want to include in a metabolic computer 
simulation. 

Now, we come to the explicit form of the proposed rate equation. Consider a 
chemical reaction with m reactants, Ai, A2, ■ ■ • , Am and n products. Pi, P2, ■ ■ ■ , Pn- 
The reaction can be expressed as 

Ai + A2 + --- Am ^^^Vs Pl+P2 + ---Pn- (2.1) 

Here we allow the case where some of the reactants maybe identical, i.e. Ai = Aj for 
some i and j and similarly for the products. The reaction is assumed to be reversible, 
with the forward reaction velocity Vp and the backward reaction velocity Vb- We 
propose a generic equation of the form 



fl{VF,VB)Y[(l+^-§-)+MVp,VB)f[ll' 



(2.2) 



1=1 ^ ' j=i \ 3 



where [Ai] with i,...,m and [Pj] with j,...,n are the concentrations of the m 
substrates and n products respectively. Here Vp and Vb are the maximal forward 
and backward reaction velocities, respectively. Ki and K'^ can be viewed as apparent 
Michaelis-Menten parameter for each reactant, defining how close each reactant to 
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saturation. The two related functions /i and /2 are introduced to take care of 
production inhibition in an exphcit way. They have the following three properties: 



The first property is a normalization condition so that at very low concentrations 
of A's and P's, the denominator is unity. The second property express the fact 
that if the reaction is only backward (i.e. Vp is zero so that the first term in the 
numerator vanishes), the rate should not be affected by the concentration of A. 
Therefore, the first term in the denominator has to be zero also. The third property 
has the same requirement as the second property but applies for the case where 
the reaction is irreversible in the forward direction (i.e. Vb is zero so that the 
second term in the numerator vanishes). In short, there is no immediate product 
inhibition directly from the enzymatic reactions for irreversible reactions where 
the product does not influence Vp and Ki. It is possible that a reaction can be 
irreversible and yet also be inhibited by its product. However, for the reactions 
that we discussed here, we are concern with irreversibility in the last elementary 
reaction where the enzyme product complex disassociates to give the free enzyme 
and product. We shall show clearly how irreversibility in our proposed generic rate 
equation and corresponding lack of product inhibition matches exactly the case in 
the full enzymatic rate equation when discussing the three cases below. We can 
also model the behavior of the enzyme in situations where product inhibition or 
activiation occurs in an irreversible reaction by specifically including a product 
modifier term. This is discussed in Sec. |5.4[ Other properties of the proposed rate 
equation will be discussed in the rest of paper. 

We emphasize that /i, /2, Vj?, Vg, {Ki}, and {Kj} should be understood as 
functions of concentrations {[^i]} and {[J'j]}, and other parameters related to their 
conditions in a milieu. This is important when we show that the generic rate equa- 
tion can be derived from the mechanistic one in Sec. [3] However, since the aim is 
to arrive at a generic rate equation that is convenient and suitable for modelling, in 
practical situations, we can assign constant values to /i, /2, Fp, Fb, {K^}, and {K'j} 
under certain physiological conditions. The number of kinetic parameters needed to 
model the system is greatly reduced. This is discussed further in Sec. |5.5[ 

Mathematically it may be ended here: The proposed rate equation has already 
been explicitly written down. It can now be used by a modeler in metabolic kinetic 
modelling. Nevertheless, the biological imphcations should be further explored to 
give a better understanding of the proposed generic rate equation. In the following 



fl{VF,VB) + f2{VF,VB) 
fliVF = 0,VB) 
/2(V>,Fs -0) 



. 







1 



(2.3) 
(2.4) 
(2.5) 
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we will first discuss its connections to three important enzymatic reactions [32) : 

Vp 

(1) isomerization: A ^ P 

Vb 
Vf 

(2) binary reaction: A + B ^ P 

Vb 

Vf 

(3) condensation-splitting: A + B ^ P + Q 

Vb 

We will explicitly demonstrate that in those situations once the rigorous kinetic 
equations are known our proposed equation can exactly reproduce them. Then we 
discuss the connections to various specific examples. 

3. Rigorous Derivation of Generic Rate Equation 



In this section we show that our proposed rate equation, Eq. (2.2 1, is exact by explic- 
itly discussing three important enzymatic reactions from the classical formalisms of 
Michaelis-Menten or King-Altman. 

3.1. Isomerization or Uni-Uni mechanism 

We start with the best known reaction, the reversible Michaelis-Menten mechanism 



(3.1) 



/c_i_2 

E + A ^ EA ^ E + P, (3.1) 

where /c+i, fc-i, fc+27 ^-2 are rate constants for each elementary reaction. It has the 
rate equation (see Eq. (2.41) in [33 , and the explanation of each parameter there) 
as 

NiA - N2P 



I + D1A + D2P 

where 



(3.2) 



^1 = I. 11 ; ^2 = , ' ° 3.3 

fc_i -t- k2 K_i -I- k2 

Cq represents the total enzyme concentration. To simplify the notation, we use 
A instead of [A] (and P instead of [P] ) to represent the concentration of A when 
writing rate equations. This simplification will be used throughout this section (Sec. 

i- 



Eq. (|3.2|) can be express exactly in the form of Eq. (|2.2|, with the constraint 

(3.5) 



that fi + J2 — Ij because Eq. (3.2 1 can be written as 

KaN,^^-K'pN2^^ 



/l(l+f^)+/2(l+tP)' 
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with 



Ka = 



/i 



Vf = 



Di 



K'p = 



Vb = 



(3.6) 
(3.7) 



D2 

N2{i-h) 

D2 

If we further choose /i = F|/(F| + V^) and = V'^/{V^ + Fj) (c./. Eq. ( [KS l 
), -fCyi, K'p, Vp and Vb are exphcitly expressed in terms of A^i, -Di and D2 as 

N^Di _ . _ NID2 



Ka = 



Vf = 



NlDl + NiDl 
NiN^Di 



Vb - 



N2NID2 



NlDl 



NlDl 



(3.8) 
(3.9) 



Thus, the rigorous equation of Eq. (3.2 1 is now exactly transformed into the pro- 



posed equation, Eq. (2.2 1, in the isomerization 



Vf^-Vb 



K' 



/l 1 



A 
Ka 



/2 1 



P 
K' 



(3.10) 



For this isomerization case, it appears there is not much to be gained in the 
reduction of parameters: One starts with five parameters, fc±i, fc±2, eg, and ends with 
four, Vp, Vb,Ka, K'p. We remark that the term "parameters" should be understood 
as "functions" in the present paper, if the dependence on temperature, pH value, 
etc, are considered, which will not be explicitly spelt out here. One would therefore 
wonder whether or not this would be a pointless game. We believe it is not. In 
writing the rate equation according to the proposed form, two features important 
in modelling kinetics dominated by enzymatic reactions are brought out clearly and 
explicitly, in addition to the stoichiometry encoding the molecular conservation law. 

The first feature is the separation of velocities of reactions, Vf,Vb, from the 
equilibrium or steady state concentrations determined in most cases by Michaelis- 
Menton like constants, Ka, K'p. Inside a cell, though concentrations of metabolites 
vary according to different physiological conditions, such variations are not very 
large. On the other hand, the reaction velocities can be strongly regulated, and can 
be changed over 6 orders of magnitude or more. For example, a reversible reaction 
can be effectively changed into an irreversible one. The proposed rate equation. 



Eq. (2.2) or Eq. (3.10 1, present this feature explicitly by separating the contributions 



into individual terms, Vf,Vb, A/Ka and P/K'p. 

The second feature is on irreversibility due to enzymatic effects. All the con- 
centrations are typically finite inside a cell, but the specific product inhibition in 
an irreversible reaction sometimes are not welcome. For example, if one naively 
switches off the backward reaction by setting iV2 = in Eq. (3.2 1, there is still an 



apparent strong product inhibition as the product concentration would still appear 



in the denominator of right hand side of Eq. (3.2 ). The subtlety lies in the definition 



Eq. (3.3 1. If irreversibility means the rate constant A;_2=0, then D2 is zero as well so 
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that no product inhibition directly from this enzymatic reaction occurs. This is not 



immediately apparent to a modeller who wants to use the rate equation Eq. (3.2 1 
but vary A^i, Di and D2 independently. Nevertheless, such apparent undesir- 
able feature is not in the proposed rate equation. The well known Michaelis-Menten 
equation for irreversible reactions without product inhibition is readily recovered, 
when the backward reaction velocity Vb = [T4t [34] : 



VfA 
Ka + A 



(3.11) 



Finally, the proposed form can be reduced back to the form for elementary 
reversible reaction without enzymatic effect. This happens when the saturation 
concentration is infinite, or A/Ka <C 1 and P/K'p <C 1. In this case the generic 
rate equation recovers the known form at low substrate concentration: 



V = apA — ubP , 



(3.12) 



with the reaction coefficients ap = Vf/Ka and as — Vs/Kp. This limiting case 
demonstrates that the proposed generic rate equation is versatile. 

3.2. Binary reaction or Bi-Uni random mechanism 

This is a typical reaction that occurs when two metabolites form another one with 
the help of enzyme(s). A typical example of such reaction is shown in Fig. [l] 



E + A 
+ 

B 

K 

EB + A 



EA 
+ 

B 

N 

k-3 

EAB==EP 



E + P 



Fig. 1. Bi-Uni Mechanism 



The steady state rate equation for the Random Bi-Uni mechanism (see Eq. (9.70) 
in [35]) is 



{NiAB + N2A^B + NsAB^ - N^P - N^AP - NqBP) 



(3.13) 



{Di + D2A + D3B + D4AB + D^A^ + DqB^ + DtA^B H- DfiAB"^ 
+DqP + DioAP + DiiBP + D12ABP) 

where the 18 terms constructed out of 11 parameters (10 rate constants: 
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fc±ij k±2, fc±3, fc±4, fc±5 and eo representing total enzyme concentration) are 
iVi = {kik^2k3k5 + k^ik2kik^)eo 
N2 = kik^kik^eo 
N3 = k2k-ikikzeQ 

Ni = (/s_ifc_2fc-3fc-5 + fc-lfc_2^-4^-5)eo 

iVs = k^ik^skik^^eo 

Ne = k^2k-3k-ik-r,eo 

Di = A;_iA;_2fc-3 + fc_ifc_2fc-4 + fc-ifc_2fc5 

D2 — kik^2k~3 + fcifc_2fc-4 + fcifc_2^5 + k^ik-ski + k^ik^k^ 

D3 = k_ik2k^3 + fc_ifc2fc_4 + k_ik2k^ + A;_2fc3fc-4 + k^2k3k^ 

D4 — kik^2k3 + fc_ifc2fc4 + kiksk-i + k2k-3k4 + k3kj^kQ 

-D5 = A:ifc_3fc4 + kik^k^ 

De = ^2^3^-4 + k2k3k^ 

Di = kikski 

Ds = ^12^3^4 

Dg = /c_i/j_2fc-5 + fc_ifc_4fc_5 + fc_2fc-3fc-5 
-Dio = k^ik^k^^ + fc_3fc4fc_5 
I?!! = k_2k3k-5 + ^3^-4^-5 

i?12 = A:3fc4fc-5 (3.14) 



A, S and P are concentrations in Eq. (3.13 1, and E in Fig. [T] represents the enzyme 



influence. This is a rather comphcated expression, though rigorously defined. 



To make the connection to the proposed rate equation, Eq. (2.2 1, the numerator 
term can be rewritten as 

AB {Ni + N2A + N3B)/Di -P (iV4 + Nr,A + NqB) /D^ (3.15) 

^ V ' ^ V ' 

VfKKaKb) Vb/K'p 

where the underbraces notation means we are rewriting (iVi + 7V2^ + -/V3i3)/£'i into 
Vf/{KaKb) and {Ni + N^A + Nf^B)/ Di into Vs/K'p. After dividing all the terms 
in the denominator of Eq. ( 3.13[ ) by Di, the denominator can be written as 



D2 D3 Di 



{l + A {D2 + D5A) +B {D3 + DeB) +AB {D4 + DjA + DsB) 

+PDa + l{DioAP + D11BP + D12ABP)) . (3.16) 

Next we introduce /i + /2 = 1 into all "l"s in the denominator above and grouped 
the constants into new functions (of A and B) in D2, D3 and D4 as shown in the 
overbraces: 

/i + AD2 + BD3 + ABDi + h{DwAP + DiiBP + D12ABP) 

+/2 + PDq + f2{DwAP + DiiBP + D12ABP) . (3.17) 
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Regrouping terms and defining x, y, z and K'p: 



y 



/i 1 + ^ 



Ih 

fi 



-B 



i + p 



/2 



DwA + DiiB 



-AB 



fi 



DuP 



D12AB 



(3.18) 



l/K'p 

and introducing two functions a and (3, we can rewrite the above as 

/i {l + Ax + By + ABz{a + /?) + ABz{l - a -/?))+ /a (l + P/K'p 
The purpose of this is so that we can factorize into 



(3.19) 



1 



A{x + BI3z) + B{y + Aaz) + ABz{l 
A{x + B(3z))(l + B{y + Aaz))+f2[l 

1/Ka I/Kb 

with Ka and Kb ss, defined provided 

(a; + Bl3z){y + Aaz) = z(l 



12(1 + P/K', 

P/K'p 



a 



13). 



(3.20) 



(3.21) 



Note that A, B and P are positive being concentration and all the Di are also 
positive, so x,y and z are positive numbers dependent on Di and P. As a result, 
the sign of the term Bj3z follows the sign of /3. Similarly, the sign of the term Aaz 
follows the sign of a. If xy = z, then a = (3 = Q. li xy < z, then a solution always 
exist with a > and /3 > 0. If xy > z, then we have a solution with a < and 
/3 < 0. We have thus shown that 



Vf 



A B 
Ka Kb 



Vb 



K' 



1 



P 

K'p 



(3.22) 



for the binary or Bi-Uni random mechanism. This is in the form the proposed rate 



equation, Eq. (2.2). This form is considerably simpler in appearance, in addition 



to having the explicit separation of the velocities of reactions and steady state 
concentrations, which in most cases are represented by the Michaelis-Menten like 
constants. Again, for the irreversible case, if we take fc_5 = 0, then A'4, N^, N^, 
Dg, DiQ, Dii and £'12 are zero in Eq. (3.141. Consequently, the original rate equa- 



tion Eq. (3.13) shows no product inhibition same as our proposed rate equation, 
Eq. pi2l. 



3.3. Condensation- Splitting or Bi-Bi ordered mechanism 

Next we consider another well known enzymatic reaction with two substrates and 
two products. 
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E + A 


"1 

^f=^ EA 


+ 


4 




Q 


E 

1 


\ 








^ 


k.3 _ . 





EQ + P ^ , " EAB==EPQ 

Fig. 2. Bi-Bi Ordered Mechanism 



The steady state rate equation for the random order Bi-Bi mechanism is very 
comphcated with 48 terms in the denominator. The steady state rate equation for 
the compulsory-order ternary-complex mechanism is simpler with 11 terms in the 
denominator and is illustrated in Fig. |2] As shown, the enzyme E must bind first 
to substrate A before it can react with substrate B. Product P is always released 
before product Q. The full rate equation is written (see Eq. (7.4) in [33] or Eq. (9.7) 
in [35]) as 

(^1^^ - N,PQ) 

{Di + D2A + D3B + D4AB + D5P + DeQ + DrBQ + DsAP ^ ' ' 

+DqPQ + DwABP + DnBPQ) . 

The 13 terms are constructed from 8 rate constants and include cq as the total 
enzyme concentration: 





= ^1^2^3^460 




= /c_ifc_2fc-3fc-4eo 


Di 


= A;-i(fc_2 + ^3)^4 


D2 


= ki{k-2 + kzfki 


D3 


^ k2k'iki 


Di 


= kik2{k3 + ki) 


D5 


= A:_ifc_2fc-3 


De 


= A;-i(fc_2 + ^3)^-4 


D7 


— k2kzk^i 


Ds 


= /cifc_2fc-3 


D9 


= (fc_i + fc-2)fc-3fc-4 


Dw 


= /Cifc2fc_3 


Dn 


= ^2^-3^-4 



(3.24) 



A, S, P and Q are concentrations in Eq. (3.231 
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To connect such complicated equation to the proposed equation, Eq. (2.2 1, we 



first divide again the numerator and denominator by Di and redefine the coef- 
ficients Ni, N2, D2, ■ ■ ■ , Dii as divided by Di. For the numerator, we can write 
A^i = Vf/{KaKb) and N2 = Vb/{K'pK'q). For the denominator we can introduce 
the functions /i and /2 subject to /i + /2 = 1 and define x,y, z,x' ,y' and z' as done 
in the Bi-Uni case to get 

X y z 



1 



- /2(1 



D2 

fi ~ 

/2 



^DsA 



-B{ 



h 

DwAB)+Q{ 



DjQ + DiiPQ)+AB{ 



/2 



DrB) +PQ{ 



fi ' 
D9 



/2 



DwP) 
-DnB) 



(3.25) 



y 



By introducing two functions a and j3 for the /i bracket part and two functions 
a' and f3' for the /2 bracket part, we can use the same method as in the Bi-Uni case 
to rewrite the denominator into 



(3.26) 



fi ( 1 

The rate equation is then 




Vf 



A B 
Ka Kb 



-Vb 



K' K, 



Q 



(3.27) 



precisely in the form of the proposed equation, Eq. (2.2 1. Furthermore, if the product 
releasing steps are irreversible, i.e. = and fc_4 = 0, then N2, D^, Dq, Dj, 
Dg,, Dq, DiQ and Dn are zero in Eq. (3.241. Hence there is no product inhibition 



in original rate equation, Eq. (3.23) same as our proposed form, Eq. (3.271 when 



irreversible. Now, there should be no question that the proposed equation is simpler, 
though formally it is equivalent to the rigorous equation as shown. 



4. Cooperativity for the case n P 

In this section, we discuss how to include cooperative effects [36] in enzyme. Hill's 
equation is often used in enzymatic kinetics to account for cooperativity. It came 
from early attempt by Archibald Hill to describe the sigmoidal O2 binding curve 
of haemoglobin in 1910 37J. If the binding of ligand by a protein with several sites 
can be adequately described by the following equation 

E + nA^E{A)n (4.1) 

then the equilibrium expression is given by 
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where Kd is the disassociation constant. The fraction of hgand binding sites filled 
is 

Y = [^(-^)"] . (4.3) 



Substituting in Eq. (4.2 1, we get 



Y = [^]" ^ . (4.4) 

We can define = Kq^ so that has units of concentration and when 1^0.5 = 
[A], half of the binding sites are filled. The reaction rate for the Hill equation [37] 
in the irreversible case then follows 

Vf[AY 



+ [^]" ■ 

This can be generalized to the reversible case [3H] 

h-l 



*^^Ao.5 ^Bp„ J [^Ao.s + Po.5 J 



1 



( [A] , [P] V 

\Ao.5 ^ Po.5 J 



(4.5) 



(4.6) 



where h is the Hill's coefficient. In our case, the proposed generic rate equation 
would be written as 



MVpM) ( 1 + (jg) ) + h(Vp, Vb) (i + a) 



Cooperativity is included by raising A/Ka to a higher power. If desired, we can also 
raise P/K'p to a higher power if cooperativity is present in the reverse direction. 
Using the same method of grouping terms and introducing additional functions in 
the previous section (Sec.|3|, it should be possible to rewrite our proposed equation 
to the desired cooperative form. 



5. Discussions 

5.1. Ansatz for fi and j-^ 



Under the conditions of Eqs. (2.3 2.5 1, a useful ansatz for the functions /i and /2 
would be 



A(T^P,F.)^ ^^,,;^^,. , (5.1) 

f^(y^^^^)-SWf^' (5.2) 

with positive numerical constants 5, e,7i and 72. The relative weight of forward 
and backward reactions are determined by S and e, while 71 and 72 determine the 
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relative influence of Vp and Vb to the generic rate equation, Eq. (2.2 1. The larger 



the 7's are, the smaller their contribution to the denominator when the reaction 
velocity is small. When 7's are infinite, there is no contribution from the smaller 
reaction. For example, if 71 = 72 = 00 and Vp > Vb, then /i = 1 and /2 = 0. On 
the other hand, the smaller the 7's are, the larger relatively the infiuence from the 
small reaction velocity. When 7's are zero, the denominator is independent of the 
reaction velocities: /i(Vf, Vb) = 5/(5 + e), /2(Vf, Vb) = e/{S + e). 

If we wish the relative effect of smaller reaction velocity to be on the smaller 
side, the numerical parameters in Eqs. (5.1 5.2 1 may be chosen as 71 = 72 = 2. 



We further take 5 ^ e = 1 so that forward and backward reactions are equally 
important. The reaction rate would then be generally written as 



VpU 



K, 




(5.3) 



For irreversible reactions, say Vb = 0, the effect of the products concentrations 




drops out naturally from the denominator. Eq. (5.3 1 can be generalized to arbitrary 
reactions. 



5.2. Thermodynamics 



We next discuss how the generic rate equation, Eq. (2.2 1 also satisfy thermodynamic 
constraints. The equilibrium constant is 



m 



(5.4) 



At equilibrium, the reaction rate is zero. Hence i^^q can also be expressed as 

n 



VbX{k, 



(5.5) 



(also called Haldane equation). Consequently, we can write the generic rate equa- 



tion, Eq. (2.2 I as 



Vp\{ 




MVp, Vb) n (1 + ^) + MVf, Vb) fl ( 1 + ^ 



(5.6) 
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where 

n 

r = ^ (5.7) 

i=l 

is often called the mass action ratio. It can be seen that the equation obeys ther- 
modynamics, since when T — K^q, the reaction rate is zero. When T > i^oq, the 
reaction is in the reverse direction, which restore the equilibrium. For T < K^q, the 
rate is in the forward direction. 

It should be pointed out that under living conditions, the medium in which 
reactions take place is active, that is, it is an open system. Since the reactions 
are usually very complicated, often not all reactants are written out explicitly. For 
example, some abundant reactants such as water (H-zO) and carbon dioxides (CO2) 
may be omitted from a reaction. In addition, some tightly regulated metabolites, 
such as adenosine triphosphate (ATP) may be regarded as constant through out the 
processes under consideration, and hence may be omitted from the reaction. Because 
of such active and open nature under living conditions, simple thermodynamical 
consideration may be misleading simply because the system is not in an equilibrium 
typically defined for a closed system. 



5.3. Comparison with other proposals 

The need to use simple (and even approximated |391 140j ) rate equations have been 
noted for some time |31j . An equation different from ours has been proposed for the 
bi-substrate case [5^ |3T] : 



1 - 



[P][Q] 

K,^[A][B] 



V[A][B] 
KaKb 



1 



Ml 

Ka 



M 



1 



M 

Kb 



K' 



(5.8) 



K may be viewed as apparent Michaelis constant with one product present. For 
example, could be estimated as the apparent Michaelis constant for A with B 
and P at appropriate concentrations and Q absent [31 . 



Comparing Eq. (5.8 1 with our propose Eq. (2.2), the first apparent difference is 



in the numerator. However, a close inspection of the rate equation, Eq. (5.8) indi- 
cates that the denominator may overestimates the contribution from the backward 
reaction. For example, when the backward reaction velocty is zero, we should not 
expect that the denominator will be directly affected by the product: it is irre- 
versible. The equation we proposed, Eq. (2.2 1 does not contain this problem since 
fi is zero whenever Vb is zero. 

Another point to note is that the derivation of Eq. (5.8 ) depends on the binding 
scheme used. Two different formulations can be obtained depending on whether 
substrate A and product P binds to one site and B and Q binding to the second 
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site, or whether A pairs with Q and B with P [29^. This requires some knowledge of 
the enzyme binding mechanism which diminish it's usefulness as a generic equation. 

The proposed generic rate equation in Ref. |29j was compared to the original 
mechanistic rate equation in the case of Bi-Bi ordered mechanism. The fit was found 
to be good under some conditions (see Fig.js]). In contrast, our generic rate equation 
can represent exactly the original mechanistic rate equation as illustrated in Fig. [3j 



0.7- 



BT 2 other Rate Equation 

Rigorous Rate Equation 
I I Generic Rate Equation 




Fig. 3. Comparison of different rate equations. Here we take the Bi-Bi ordered mecfianism witli 
8 rate constants as an example. In Ref. |29l . a generic bi-substrate rate equation was fitted to 
the mechanistic rate equation according to a set of data points for the substrates and products 
concentrations (Conditions CIII in 1291 . see their Fig. 2). The fit at product concentration P = 
Q = 0.1 is considered good. Nevertheless, our proposed generic rate equation, Eq. | |2.2| , is capable 
of reproducing exactly the original rigorous mechanistic rate equation. We note that it is likely 
the other equation, Eq. | |5.8| may also be cast into the rigorous rate equation using the method 
shown in Sec. [s] 
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5.4. Adding modifiers 

Modifiers are substances that interact with enzymes to either increase or decrease 
their catalytic activity. Inhibitors diminish the rate of reaction while activators in- 
crease the rate. These are very important in the regulation of metabolic pathways 
in organisms. For reversible inhibitors, they can be classified into competitive in- 
hibitors and noncompetitive inhibitors. Competitive inhibitors compete with the 
substrate for binding to the enzyme forming an enzyme-inhibitor complex. The in- 
hibitor does not act on the enzyme-substrate complex once it is formed. For the 
simple Michaelis-Menten kinetics, the resulting form is 



Vf[A] 



(5.9) 



Vp is not affected but Ka is increased. Non-competitive inhibitors bind with the 
same affinity to the enzyme and enzyme-substrate complex. It does not change Ka 
but decrease Vp- 



-Il]/Ki 



Vf[A] 

(5.10) 



Ka + [A] 

When the affinities are different, the inhibition is called mixed. The net result is to 
increase in Ka and decrease in Vp- A special case occurs if inhibitor only binds to 

the enzyme-substrate complex but not to the enzyme, this is called uncompetitive: 

(^)''H^l ,,,,, 



(i+ItW) + [-^1 ' 

To include the effects of inhibitors and activators, we can multiply the unmodi- 
fied rate equation with a function of the sigmoidal form 

^ (5.12) 



for an inhibitor, / and 



l + ([/]/a)2 

m/br 

i + m/by 



(5.13) 



for an activator, T. The parameter a can be adjusted to match the concentration of 

the inhibitor at which inhibition becomes significant. The same applies for b in the 

case of an activator. Other power in the denominator besides 2 can also be used. 

It is clear that when Vp is modified by multiplying the modifier function, this is 

similar to the case for non-competitive inhibitor in Eq. (|5.10 1. For the compet itive 

' ' 

K 



inhibitor suitable choice of 1/Ki = (1 -I- j^)/a substituted into Eq. (|5.9l 
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Vf[A] 
Ka (1 + [I]/Ki) + [A] 
Vf[A] 



KA{l + i[I]/a)il+[A]/KA)) + [A] 
Vf[A] 



{l + [I]/a){KA + A)- 



(5.14) 



This is the same form of the inhibitor that we used in Eq. (5.121. Hence we have 



shown that both competitive and non-competitive modifiers can generally be ac- 



counted for by the modifier form that we used in Eq. (5.121 



5.5. Uses of proposed rate equation 

We have demonstrated above that in situations where exact enzymatic reactions 
are known, the proposed rate equations are the same as them. In realistic modelling 
situations, the enzymatic information is typically incomplete. The proposed equa- 
tions are easy to use because of its simplicity and its parameters or functions are 
easy to interpret. We elaborate on this observation further below. 

Earlier in deriving the generic rate equation (in Sec. [s]), we take Vp, Vb and, 
i^i's where i are metabolites in the enzymatic reaction, to be functions of concen- 
trations of metabolities and rate constants in the elementary reactions. We already 
alluded to the situation of incomplete information, and to the nearly impossibility 
of obtaining all the relevant in vivo kinetic parameters experimentally. Fortunately, 
it is well known that the biological networks, either genetic or metabolic, or others, 
are typically robustly built [21 |6l |8j [HI [THJ [17] at functional level. Therefore, very 
detailed kinetic information may not be necessary in practical situations where we 
want to simulate a metabolic model under some physiological or functional condi- 
tions. We may assign typical values, which may be approximated by constants, to 
the parameters Vp, Vb and Ki in the simulation. Since the operational meaning of 



these parameters are clear (see below Eq. (2.2)), we can adjust these parameters in 
a meaningful way to match physiological conditions. 

There is another way to use the proposed enzymatic equation: The compact 



form of Eq. (2.2 1 with approximated values in the kinetic parameters may serve as a 
starting point to introduce empirical description. Once the values of parameters are 
fixed, specific experiments to measure them can be performed, and the comparison 
between theoretical and experimental efforts can lead to better understanding and 
further study. Clearly, this appears a more straightforward procedure than tuning 



the rate constants of the elementary reactions (such as the fc^'s in Eq. (3.14 1) in the 



original rate equation (such as Eq. (3.13 1). Effectively, the number of parameters in 
the model is reduced. Can such a strategy be useful? Our experience so far suggests 
yes. 

From a different angle, a third possible use of the proposed enzymatic equation 
is to regard it as a format to document kinetic parameters from experiments: The 
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form is formally exact, with easy to interpret parameters, and, is ready for compu- 
tational modeling. In particular, it would be of great value if experimental results 
on reversible reactions, particularly on the backward velocities, are documented. 



We have been using the generic rate equation, Eq. (5.3 1 to construct a kinetic 
model of the metabolism of Methylobacterium extorquens AMI |4Tl|42l|43] consist- 
ing of about 80 reactions and 80 metabolites. The number of metabolites are the 
minimum in order to understand the metabolic kinetics at a global level [41 j l44] . We 
are able to find a set of parameters {Vp, Vb and Ki's) that give a steady state solu- 
tion |45j . Our results indicate that the model is robust with respect to variations in 
kinetic parameters. This is in agreement with observations that quite large errors in 
kinetic parameters for most of the enzymes have very little effect on the calculated 
flux [21]. This result is very encouraging as it suggest even larger metabolic network 
can be built on the same structure for kinetic equations. 

Reducing the number of kinetic parameters are also useful when metabolomic 
data are available for analysis. Under such situation, one could attempt to fit the 
kinetic parameters to the metabolomic data. One way may be to use Genetic Algo- 
rithms to optimize the parameters |46j . Simplifying the number of kinetic parame- 
ters cuts the computational effort and makes the problem much more tractable. 



6. Conclusions 

We have proposed a generic kinetic rate equation. It is shown to be formally identical 
to rigorous equations and has a considerably simple and generic form. It is already 
in a form workable for arbitrary number of substrates and products with different 
stiochiometry. Modifier functions can be added to reflect inhibition and activation. 
Cooperativity can also be included. 

In order to do kinetic modelling according to rigorous equations, numerous pa- 
rameters need to be measured. Instead, the proposed generic rate equation reduces 
the required number to the most essential ones, both forward and backward reac- 
tion velocities, Vp, Vb, and the Michaelis-Menten like constants, K,K', for both 
substrates and products, which have clear chemical and biological meanings. Even 
if those essential parameters are not always available under appropriate physiolog- 
ical conditions related to living organisms in databases, their adjustments during 
modelling are much more managable. 

The kinetic parameters from enzyme characterization is also dependent on the 
rate equation being tested. There is an effort underway to ensure uniformity in 
reporting functional data [471 . The proposed rate equation with only essential pa- 
rameters may aid such effort to have enzyme characterized. With a database of 



enzyme parameters for a generic rate equation of the form Eq. (2.2 1, the result- 
ing uniformity will make reconstruction of genome scale metabolic networks using 
kinetic models will become much more feasible. 
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